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• Using Computational Fluid Dynamics (CFD) to 
predict a flow field is an approximation to the 
exact problem and uncertainties exist. 




• There is a method to approximate the errors in 

CFD via Richardson's Extrapolation. 

— This method is based off of progressive grid 
refinement. 


• Unless using a Structured Grid with every other 
point, some interpolation method must be used. 
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Summary of Richardson's Extrapolation jc~r> 
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• Navier Stokes Equations 

— 2 nd order, non-homogeneous, non-linear partial 
differential equations 

• Richardson's Extrapolation is used to produce 
4 th order accurate solution from separate 2 nd 
order accurate Navier Stokes Solutions 
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Summary of Richardson's 
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• ASME V&V 20-2009 Outlines a 5-step Procedure 
to Richardson's Extrapolation using Roache's 
(1998) Grid Convergence Index (GCI) Method 

• Assumptions 

1. Three discrete solutions are in the asymptotic range 

2. Meshes have a uniform spacing over the domain 

3. Meshes are related through systematic refinement 

4. Solutions are smooth 

5. Other sources of numerical error are small 
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ASME V&V 20-2009 5-step Procedure to 
Richardson's Extrapolation 
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• 5 Step Procedure for Uncertainty Estimation 

- Step 1: Representative Grid Size 


// 





where 

N = total number of cells used for the computations 
A V. = volume of the i th cell [4] 

h x < h 2 < h 3 
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ASME V&V 20-2009 5-step Procedure to 
Richardson's Extrapolation 
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- Step 2: Select 3 significantly (r>1.3) different grid 
sizes , /7 

= hJh 


21 


32 


V 1 


— Use CFD Simulation to analyze key variables, Sr 

*32 = V3- ?! 

*21 = <P 2 - <Pl 



JOHN F KENNEDY SPACE CENTER 


ASME V&V 20-2009 5-step Procedure to 
Richardson's Extrapolation 



- Step 3: Calculate observed order, p 


p = [1 /ln(r 2l )][ln\s 32 / e 21 \ + q(p )\ 


q (p) = In 


s = 1 • sign(£ 32 / s 21 ) 


(H’i - 

“ S \ 

V 32 - 

-sl 
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ASME V&V 20-2009 5-step Procedure to 
Richardson's Extrapolation 
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- Step 4: Calculate extrapolated values 
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ASME V&V 20-2009 5-step Procedure to 
Richardson's Extrapolation 
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Step 5: Calculate Fine Grid Convergence Index & 
Numerical Uncertainty 


GCI f t 


The Factor of Safety, Fs = 1.25 


Fs • e 


21 


a 


yP 

' 21 


1 


- Assumption that the distribution is Gaussian about the fine grid, 90% 
Confidence 


u 


num 


= GCI / 1 .65 


Solver Interpolation 
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• FLUENT 

• Includes a Mesh-to-Mesh Interpolation 

• Performs a zeroth-order (nearest neighbor) interpolation 

• Designed for initial conditions from a previous solution 




• OPENFOAM 

• Mapfields fuction interpolation 

• Used for initialization of a solution from a previous model 


• Using these 'zeroth-order' interpolation schemes is not 
sufficient for comparing errors from the mesh 
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Matlab Interpolation Schemes 



• Matlab 

• High level language used for numerical computations 

• CFD data is in various forms 

• ID, 2D, 3D, uniform, non-uniform 

• Generic Scheme is sought for all CFD data 

Matlabifunction 

interpl interp2 interp3 

InterpolationH/lethod 
'nearest'HQNearestfflieighborSInterpolation 
'linear'SHinearinterpolation[|default) 

'spline'SSCubicBplineQInterpolation 
'pchip'HIFiecewiseilubiclHermiteQInterpolation 
'cubic'0 

'v5cubic'0Hubic0nterpolationlSjsedl3nll/latlabE& 
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X 

Xl|uniformly-spaced0)nly) 

XI|uniformly-spacedl2bnly) 

X 
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Example Problem 



- Fully developed flow between parallel plates 

• Exact Solution to Navier Stokes 


• Provide a good example of errors that can be induced 
from interpolation 
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al|m)0 

0.10 

rhol| kg/m 3)0 

1.2250 

mul|Ns/m2)0 

0.000017890 

dp/dx|N/m3)0 

-0.0040 
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Example Problem 
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- Constructed a 
CFD Model in 
FLUENT 

• 3 Grids 

- Coarse, 7,140 
Cells 

- Medium, 14,186 
Cells 

- Fine, 24,780 
Cells 


CFD vs. Exact 



Velocity (m/s) 


^EXACT 
• Coarse 
▲ Medium 
■ Fine 
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Example Problem 
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• Interpolation Direction? £ v 



LAUNCH SERVICES 


<e 2 ~ <p, 

«f> 3 - <P 2 


1. Interpolate Coarse and Medium Mesh -> Fine 



1. Interpolate Medium and Fine Mesh -> Coarse 
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Position Y (m) 


Example Problem 



JOHN F KENNEDY SPACE CENTER 


1. Linearly Interpolate Coarse and Medium Mesh -> Fine 
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<P 2 

medium 


Extrapolated Solution on Fine Grid 

0.0600 



0.0000 0.0050 0.0100 0.0150 0.0200 0.0250 0.0300 
Velocity (m/s) 


■ Extrapolated Solution 
-•-EXACT 
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Max % Error 
Extrapolated 
Values 

Average % 
Error 

Extrapolated 

Values 

0.8950 

0.0596 
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Position Y (m) 


Example Problem 
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2. Linearly Interpolate Fine and Medium Mesh -> Coarse 
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Extrapolated Solution on Coarse Grid 


0.0600 

0.0500 

0.0400 


/ 



Velocity (m/s) 



UCF 


<p. 


3 * 

coarse 



Max % Error 

Average % 


Extrapolated 

Error 

■ Extrapolated Solution 
“—EXACT 

Values 

Extrapolated 


Values 


0.0792 

0.0175 
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Example Problem 


Interpolation Direction? 



1. Interpolate Coarse and Medium Mesh -> Fine 


Extrapolated Solution on Fine Grid 



■ (itr«poUtrd Solution 
“♦"EXACT 


Velocity (m/%) 


Max % Error 
Extrapolated 
Values 

Average % Error 
Extrapolated 
Values 

0.8950 

0.0596 


11 Interpolate Medium and Fine Mesh -> Coarse 


Extrapolated Solution on Coarse Grid 



■ t Ktr»poljtPd Solution 
“♦“EXACT 


Velocity (m/%) 


Max % Error 
Extrapolated 
Values 

Average % Error 
Extrapolated 
Values 

0.0792 

0.0175 









Example Problem 
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- Interpolating to the coarse grid was selected 

- Other interpolation methods 

• "nearest" - Fluent's Mesh-to-Mesh 


• "linear" - Matlab 

yfi = interplJfinefi^J^inei^lJ^oarsefi^J/linear') 

• "cubic" - Matlab 

yfi = interpl(fine(:,2),fine(:,l),coarse(:,2)/cubic') 
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Example Problem 



• "nearest" 



- Fluent's Mesh-to-Mesh 
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Example Problem 
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• "linear" - Matlab 

yfi = interplffinel^^fine^l^coarsefi^/linear') 
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Example Problem 
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• "cubic" - Matlab 

yfi = interpl(fine(:, 2), fine(:,l), coarse):, 2)/cubic') 
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Matlab Interpolation Schemes 



• Extending the Interpolation Schemes to 2D and 3D 

• Interp2 and Interp3 Matlab Functions 

• Require use of MeshGrid 

• Transforms the domain of vectors into arrays 

• For Meshes in the 4 million to 8 million Cell Range 

• Error "Maximum variable size allowed by program is 
exceeded" 

• Griddata Function 

• Nearest, Linear, Natural, Cubic, and v4 

• Nearest, Linear, and Natural are the only options available in 2D and 
3D 

• The only options available for ID, 2D, and 3D 

• Interpl and Griddata - 'nearest' 

#>UCF 
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• Airflow around encapsulated spacecraft 

- Matlab griddata 'linear' option used 

- Interpolating Fine and Medium Grid onto Coarse 
Grid 
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imf 3D Example 
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• Comparing using a Line Plot 
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3D Example 



• Comparing using a Line Plot 
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Conclusion / Recommendation 



• By comparing the interpolation schemes in one, two, and 
three dimensions and investigating the options that are 
readily available in Matlab 

• Recommend the "linear" option be used when comparing the 
error or uncertainty due to the grid 


• If coarse grid has the level of detail required 

• Recommend interpolating from the fine and medium grids onto 
the coarse grid 

• Lower Error in the Extrapolated Solution 

• Smaller Data Set 

— Future Work include higher order interpolation schemes 
in 3D (Radial Basis Function Interpolation, 4 th order) 


interpl or griddata Matlab commands 
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